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Abstract 

The generalized sequence of numbers is defined by Wn = pWn-i + 
qWn-2 with initial conditions Wo ~ a and Wi = b for a,b,p,q £ Z 
and n > 2, respectively. Let W„ — circ{Wi,W2, ■ ■ ■ ,Wn)- The aim of 
this paper is to establish some useful formulas for the determinants and 
inverses of W,i using the nice properties of the number sequences. Matrix 
decompositions are derived for W„ in order to obtain the results. 



1 Introduction 



The n X n circulant matrix C„ = circ(co, ci, . . . , c„_i), associated with the 
numbers cq, ci, . . . , c„_i, is defined by 



Co Cl 
c,i-i Co 

Cl C2 



Cn-2 C„_i 

Cri-3 C„_2 

Co Cl 
Cti~1 Co 



(1) 



Circulant matrices have a wide range of applications, for examples in sig- 
nal processing, coding theory, image processing, digital image disposal, self- 
regress design and so on. Numerical solutions of the certain types of elliptic 
and parabolic partial differential equations with periodic boundary conditions 
often involve linear systems associated with circulant matrices [9-11]. 
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The eigenvalues and eigenvectors of C„ are well-known [4,14]: 

n-l 

Xj = ^CkUJ^'', j = 0, l,...,n- 1, 
fe=o 

where cj := exp(2^) and i := \/^ and the corresponding eigenvectors 

Xj = . j = 0, 1. 

Thus we have the determinants and inverses of nonsingular circulant matrices 
[1,3,4,14]: 

n—1 n — 1 

det(C„) = [](^Cfe^^'^-) 

and 

C~'^ = circ{ao,ai, . . . ,a„_i) 

where Ofc = ^ S"=o and = 0, 1, ... ,n — 1) [4]. When n is getting 

large, the above the determinant and inverse formulas are not very handy 
to use. If there is some structure among cq, Ci, . . . , c„_i, we may be able to 
get more explicit forms of the eigenvalues, determinants and inverses of C„. 
Recently, studies on the circulant matrices involving interesting number se- 
quences appeared. In [1] the determinants and inverses of the circulant matrices 
An = circ{Fi, F2, . . . , Fn) and -B„ = circ{Li, L2, . . . , Ln) are derived where F„ 
and Ln are nth Fibonacci and Lucas numbers, respectively. In [2] the r-circulant 
matrix is defined and its norms was computed. The norms of Toeplitz matri- 
ces involving Fibonacci and Lucas numbers are obtained [5]. Miladinovic and 
Stanimirovic [6] gave an explicit formula of the Moore-Penrose inverse of singu- 
lar generalized Fibonacci matrix. Lee and et al. found the factorizations and 
eigenvalues of Fibonacci and symmetric Fibonacci matrices [7] . 

The generalized sequence {Wn{a,b;p,q)} or {Wn} of numbers is defined 
by Wn = pWn-i — qWn-1 with initial conditions Wq = a and W\ = b for 
a,b,p,q € Z and n > 2, respectively [16, pp. 161]. Let a and (3 be the roots of 
x"^ — px + q = 0. Then the Binet formula of the sequence {Wn} is 

Wn = 3 — 16, pp. 161 

a — p 

where A = b — a(3 and B = aa — b, a + (3 = p, afi = q and a — fi = sJ'fP' — 4(7. 
Let p = q = a = ^,b = \ and a = 2, 6 = 1, then G„ are Fn nth Fibonacci 
and L„ nth Lucas numbers, respectively. While p = 2, g = 1, if a = 0, & = 1 
and a — b — 2, then Gn are P„ nth Pell and Qn nth Pell-Lucas numbers, 
respectively. 

Let W = circ{Wi, W2, . . . , Wn)- The aim of this paper is to establish some 
useful formulas for the determinants and inverses of W using the nice properties 
of the number sequences. Matrix decompositions are derived for W in order to 
obtain the results. 
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2 Determinants of circulant matrices with the 
number sequence 

Recall that W„ = circ{Wi, W2, ■ ■ ■ , Wn), i.e. where Wk is fcth element of se- 
quence {W„}, with the recurence relations Wk = pWk-i + qWk-2, the initial 
conditions Wq = a, Wi = 6 (fc > 2) when q is positive real number. Let a and 
(3 be the roots oi — px + q — while — Aq 0. Using the Binet formula 
[16, pp. 161] for the sequence {W„}, one has 



a- P 



[16, pp. 161] 



(2) 



Theorem 1 Let n > 3. Then 



detf 



{b'' -W2Wn){b- 



n-l 



Y,[{hWk+l-W2Wk) X 



fe=2 



x{b-Wn+lt-^{qWn 



xn — ki 



where Wk is kth the element of the sequence {Wn} 

r3 



(3) 



Proof. Obviously, det(W3) = b^ + W^ + W| - 3&W2VF3. It satisfies Q. For 



n > 3, we select the matrices K and L so that when we multiply W„ with K 
on the left and L on the right we obtain a special Hessenberg matrix that have 
nonzero entries only on first two rows, main diagonal and subdiagonal: 



K := 



1 

W2 

Wi 

-q 




-p 
-q 



... 


' 


... 


1 


. . . 1 


-p 


. . . -p 


-q 


... 





... 






(4) 



and 



L := 








q(Wr 


-Wo) 


Wi- 




q{Wr 


-Wo) 


Wi- 




q(Wr 


-Wo) 



( 

( q(Wr.-Wo) \ 
\ Wi-W„+i ) 



n — 3 



n— 4 



q(W^-Wo) 
Wi-W„+i 

1 





1 

1 
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Notice that wc obtain the fonowing equivalence: 



9n W^-'-^^^ P^„_i 
Wl - Wn+1 



Wi 



n— 3 



W3 



W2 



W£ 
Wi 



Wi - Wn+l 



and M is Hessenberg matrix, where 



k=2 



\Wl- Wn+1 



n — k 



and 



5n := W, + 



^ fc=2 \ ^ ^ 



Then we have 

det(M) = det(if) det(W„) det(L) 

Since 



Wn+lY'-'gn. 



for ah n > 3, 
and follows. 



det(i^) = dct(i) = 

det(A')det(L) = 1 



1 , n — 1 or 2 mod 4 
1, n = or 3 mod 4 



3 Inverses of Wr 
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Let Ck,n be an n X n A;-circulant matrix. Then 

Ac.,„(A) = A"-'"n(^"^-%) 
i=o 

is the characteristic polynomial of the matrix Ck^n where 
and 

7% 

y= - (See [15], pp. 3). 

Since GCD{l,n) = 1, then n = m and the characterictic polynomial of the 
matrix Ci,„ is 

Ac.,„(A) = nV-%). 

Therefore, is not the eigenvalue of Ci_„. Then dct(Ci^„) 7^ for n > 3. 

We will use the well-known fact that the inverse of a nonsingular circulant 
matrix is also circulant [14, p.84] [12, p.33] [4, p.90-91]. 

Lemma 2 Let A = {aij) be an (n — 2) x (n — 2) matrix defined by 



Then A ^ = (a^^ ) is given by 



Wi-Wn+l, 1=3 
q{Wo-Wn), 1=3 + 1 
0, otherwise. 



(g(Wo-ty„))'-^' . ^ . 
(Wi-W„+i)'-i+i ' ' — J 

0, otherwise. 



Proof. Let B := (bij) = BB ^. Clearly, 6^ = Y^^^^ikO-kj- When i = j, we 
have 

If i > j, then 



similar for i < j. Thus, BB^^ = In-2- ■ 
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Theorem 3 Let the matrix W„ he W„ = circ{Wi,W2, W„) (n>3). Then 

the inverse of the matrix W„ is 

= circ{'Wi,W2, Wn) 



where 
Wi = 



W2 = 



W3 



W4 = 



Wr, 



Qn Qn 



fc=i 



J2 Wn-k - 



E 



\Wi-Wn+l 
fq{Wo-Wn 



l+P 



Wi - Wn+1 



W2 1_ 

gnWi gn{Wi - Wn+l) ^^\Wr - W^+l 

W2Wn-k 
Wl 



X Wn-k+1 - 



WiWs - Wj 
9nWi{Wi-Wn+l 
1 



gniWi 


- Wn+l) 




Q 


gn{Wi 


- Wn+l) 




1 



W4 



W2W3 

Wl 



Wn+2~W2 \ 
WiJ \ Wi-Wn+lJ 



Wi \ fiWo- w„)(w„+2 - W2) 



{Wl - Wn+lf 



- 1 



gniWi - Wn+l) 
'Wn+2-W2 



Wl - Wn+1 

Wo-Wn 
Wl - Wn+1 



Wn- 
n-2 



W2Wn- 
Wi 



W2Wn-2 , 

^"-^ " -^w~ ' 



fe=2 



Wl 



{Wo-Wn){Wn+2-W2) 
{Wl - Wn+l)^ 



n-1 



for gn:=W^ 



Wl 



+ 



E(K-^) {^^y~')- 



fc=2 
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Proof. Let 



1 - 



U 



an 

Wi 

1 



gn 



Ul3 
WiWn 



Ml4 

W2Wr, 



Ul5 

W2Wn 



W4 _ W2W3 
Qn Qn Wi 



W3 W, 
9r 







1 

1 

where 
and 



9n = T.^k(- 



fe=2 

and 



Wi - Wn+l 



Let = (iia(7(Wi, (?„). Then we can write 

KWnLU = H®A 

where H ®Ais the direct sum of the matrices H and A. Let T = LU. Then we 
have 

W-^ =T{H-'^ ®A-^)K. 

Since the matrix W„ is circulant, its inverse is circulant from Lemma 1.1 [1, 
p.9791]. Let 

:= circ{wi,W2, ■ ■ ■,Wn)- 
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Since the last row of the matrix T is 

9n 9nWi ' Qn QnWi ' Qn QnWi 

Wi W2W3 W3 W| 



9n 9nWx Qn QnWi 

the last row components of the matrix are 
W2 1 



W2 



W3 
W5 



W2Wn-k 



WiWs - w| 

1 

gn{Wx-Wn+l) 

q 

9n{Wx - Wn+l) 



Wi 



Wi\ (Wn+2-W2 



Wl J \ Wl-Wn+l 



^ _ ( {WQ-Wn){Wn+2-W2) _ \ 

' Wj\ {Wl-Wn+lY ) 



Wn = Wr, 



Wl = 



9n{Wi-Wn+l) 



W2W, 



n-1 



Wl 



f Wn+2-W2 
\Wi-Wn+l 



n-2 



fc=2 

{Wo - Wn)iWn+2 - W2) 
{Wl - Wn+l)^ 
1 



W2Wn- 



Wn-1 — I X 



fc-1 



9n 9n 

{Wl - Wn+l) 

n-2 



p[Wn 



W2Wn 
Wl 



Wl 



+ 



Wl 

Wl - Wn+l 



kl,,, W2Wn-k-l \ f Wo-Wn 

2^ g'= ( Wn-k 



fc=l 



Wl 



\Wi-Wn+l 



X 



Wo-Wn 
Wl - Wn+l 
n-1 

where gn = W,- ^ 



Wl 



is circulant matrix, the proof is completed. 



E (K^ - ^) («^)"" ) • S-e W- 

k=2 ^ ' 
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